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ABSTRACT 

This paper considers the distribution of dust which originates in the break-up of planetesimals 
that are trapped in resonance with a planet. The effect of radiation pressure on small dust grains 
causes their orbits and so their spatial distribution to be significantly different to that of the parent 
planetesimals which previous work has shown to be clumpy. It is shown that there are three distinct 
grain populations: (I) large grains (> a few mm) have the same clumpy resonant distribution as the 
planetesimals; (II) moderate sized grains (a few fim to a few mm) are no longer in resonance and 
have an axisymmetric distribution; (III) small grains (< a few fim) are blown out of the system 
by radiation pressure immediately on creation and so have a density distribution which falls off as 
T (X l/r, however the structure of these grains can be further divided into two subclasses: (Ilia) 
grains produced in the destruction of population I grains that exhibit trailing spiral structure which 
emanates from the resonant clumps; and (Illb) grains produced from population II grains that have 
an axisymmetric distribution. Since observations in different wavebands are sensitive to different sized 
dust grains, multi- wavelength imaging of debris disks provides a valuable observational test of models 
which explain structure seen in sub- mm observations as due to resonant trapping of planetesimals. For 
example, a disk with a coUisional cascade dust size distribution with no blow-out grains would appear 
clumpy in the sub- mm (which samples population I grains) , and smooth at mid- to far-IR wavelengths 
(which sample population II grains). The wavelength of transition from clumpy to smooth structure is 
indicative of the mass of the perturbing planet. The size distribution of Vega's disk is modeled in the 
light of the recent Spitzer observations showing that collisions with the large quantities of population 
III grains seen in the mid- to far-IR may be responsible for the low levels of population II grains in 
this system. The origin of these population III grains must be in the destruction of the grains seen in 
the sub-mm images, and so at high resolution and sensitivity the far-IR and mid-IR structure of the 
Vega disk is predicted to include spiral structure emanating from the sub-mm clumps. Such structure 
could be detected with MIRI on the JWST and, if so, would confirm the presence of a planet at 65 
AU in the Vega disk as well as determine the direction of its orbit. 



Subject headings: celestial mechanics - 
individual (Vega) 



- circumstellar matter — planetary systems: formation — stars: 



X 



1. INTRODUCTION 

Debris disks are the dust disks that are fed by the 
collisional grinding down of extrasolar planetesimal belts 
that are analogous to the Kuiper Belt in the Solar System 
(Backman & Paresce 1993; Wyatt et al. 2003). One of 
the most defining features of the debris disks that have 
been imaged is that most of these disks are not smooth, 
but clumpy. Such clumps are seen in sub-mm and mm 
images of the dust disks around Vega (Holland et al. 
1998; Koerner, Sargent & Ostroff 2001; Wilner et al. 
2002), £ Eridani (Greaves et al. 1998; Greaves et al. 
2005) and Fomalhaut (Holland et al. 2003), as well as 
in optical and mid-IR images of (3 Pictoris (Kalas et al. 
2000; Telesco et al. 2005). 

The origin of these clumps has been widely debated. 
One possible interpretation is that the clumps are in 
fact unrelated, and are the chance superposition of back- 
ground objects. However, this has been ruled out on 
statistical grounds based on number counts of such ob- 
jects (e.g., Holland et al. 1998; Telesco et al. 2005). 
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Another possibility that has been considered is that the 
clumps were created relatively recently in collisions be- 
tween large planetesimals; features in the structure of the 
zodiacal cloud have been interpreted in this way (Der- 
mott et al. 2002; Nesvorny et al. 2003). However, the 
high mass of material seen in the clumps in disks im- 
aged at long wavelengths indicates that such an event 
would have to have involved two planetesimals at least 
the size of Pluto, which given the expected frequency of 
such events again rules out this interpretation on sta- 
tistical grounds (Wyatt & Dent 2002; see also Kenyon 
& Bromley 2005). This mechanism remains a possibility 
for clumps seen at shorter wavelengths, for which smaller 
colliding planetesimals are required to create detectable 
clumps, for clumps seen at small orbital radii, and for 
those seen toward relatively young systems like (3 Pic- 
toris (Telesco et al. 2005; Kenyon & Bromley 2005). 

The favoured interpretation for the clumps seen in sub- 
mm and mm images is that they are associated with 
material that is trapped in resonance with a planet or- 
biting in the disk. Two classes of model invoke plan- 
etary resonances to explain the clumps, and they dif- 
fer in the mechanism by which material ends up in the 
resonances. In one model, dust migrates inward due to 
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Poynting- Robertson (P-R) drag and upon encountering a 
resonance, resonant forces temporarily halt the migration 
causing a concentration of dust in the resonances. This 
is the mechanism by which the Earth's clumpy resonant 
ring is thought to have formed (Dermott et al. 1994). 
The clumpy structures of the disks of Vega, e Eridani 
and Fomalhaut, have all been modeled in this way indi- 
cating the presence of planets more massive than Saturn 
oribitng at several tens of AU from these stars (Ozernoy 
et al. 2000; Wilner et al. 2002; Quillen & Thorndike 
2002; Deller & Maddison 2005). However, these models 
suffer from the problem that, unlike dust in the zodiacal 
cloud, inward migration due to P-R drag is not signif- 
icant in these systems because their disks are so dense 
that collisions occur on much shorter timescales. This 
means planetesimals are ground into dust fine enough to 
be removed by radiation pressure before P-R drag has a 
chance to act (Wyatt 2005). This is not the case for dust 
coming from the asteroid belt which is much less dense 
meaning that collisions are much less frequent. 

In the other model, the parent planetesimals of the 
dust were trapped in resonance with a planet which mi- 
grated outward early in the history of the system, in 
the same way Kuiper Belt objects were trapped in res- 
onance with Neptune when it migrated outward (e.g. 
Malhotra 1995; Levison & MorbideUi 2003). Such mi- 
gration could have been caused by angular momentum 
exchange when the planetary system scattered the resid- 
ual planetsimal disk (Fernandez & Ip 1984; Malhotra 
1993; Hahn & Malhotra 1999; Ida et al. 2000). Wyatt 
(2003; W03) used numerical simulations to model the 
dynamical and spatial structure of a planetesimal disk 
resulting from the outward migration of a planet. W03 
also showed how comparing this structure with that seen 
in sub-mm images of the Vega disk not only explained the 
observed structure without having to invoke P-R drag, 
but also allowed constraints to be set on the planet caus- 
ing the structure, since only a specific range of planet 
mass and migration rate can cause the observed struc- 
ture. A Neptune-mass planet which migrated from 40-65 
AU over 56 Myr was proposed which resulted in plan- 
etesimals being captured predominantly into the planet's 
2:l(u) and 3:2 resonances. One of the limitations of the 
W03 model though was that the dust seen in the sub-mm 
images was assumed to have the same distribution as the 
planetesimals. Small dust grains generally have different 
orbital parameters to their parent planetesimals because 
radiation pressure causes them to effectively see a lower 
mass star than did their parents. This means that the 
dust distribution may be significantly different to that of 
the planetesimals. 

This paper considers how the distribution of dust 
arising from a population of planetesimals previously 
trapped into resonance by a migrating planet differs from 
that of the planetesimals themselves. The distributions 
of two different types of grains are treated separately: in 
^ numerical simulations are used to consider the dis- 
tribution of grains that remain gravitationally bound to 
the star; []21models the collisional evolution of the disk to 
determine the distribution of grains that are put on hy- 
perbolic orbits by radiation pressure as soon as they are 
created. The findings are summarised in which shows 
how the disk can be divided into three grain populations 
sorted by grain size, each of which exhibiting a distinct 



structure. This section also considers which populations 
will dominate observations in different wavebands and 
discusses the implications for interpretation of the Vega 
disk in the light of the recent Spitzer images of its struc- 
ture (Su et al. 2005). The conclusions are given in |21 

2. DISTRIBUTION OF BOUND GRAINS 

Small dust grains are acted on by a radiation force 
that can be parametrised by a factor (3 which is the ratio 
of the radiation force to that of stellar gravity (Burns, 
Lamy, & Soter 1979). The factor /3 is a function of par- 
ticle diameter with smaller particles having larger values 
of (3 which for large grains falls off oc 1/D. There are 
two components to the radiation force: radiation pres- 
sure, which is the radial component that means that the 
particle effectively sees a smaller mass star by a factor 
1 — /3; and Poynting-Robertson (P-R) drag, which is the 
tangential component that makes the particle's orbit spi- 
ral in toward the star at a rate -13/a. While it 
is included in the following numerical simulations, the 
P-R drag force is not dominant in causing dust grains to 
have different distributions to their parent planetesimals. 
Rather that change in structure is caused by the fact 
that on creation the dust grains have the same positions 
as their parents, and similar velocities, but see a smaller 
mass star. This causes the parameters describing their 
orbits to be different to those of the parent planetesimal. 
Dust grains created in the break-up of a planetesimal 
that had /? = 0, and for which orbital elements at the 
time of the collision were a, e, /, fi, a), and /, move in the 
same orbital plane as the parent, 1=1 and = ^l, but 
on orbits with semimajor axes, a , eccentricities, e , and 
pericenter orientations, uj , that are given by (Burns et 
al. 1979; Wyatt et al. 1999): 

a=ail~l3)/ [1 - 2/3(1 + e cos /)/(l-e2)] , (1) 
e' = (1 - fSyWe^ + 2/3e cos f + P^, (2) 
uj — Lj = f — f = arctan [/3sin//(/3cos/ -I- e)]. (3) 

Two types of dust grains are created: those (large 
grains) with /3 < 0.5 that remain on bound orbits; and 
those (small grains) with /3 > 0.5 that are blown out 
of the system on hyperbolic orbits as soon as they are 
created. For the former, the important question to ask 
is how a particle's new orbital elements affect its subse- 
quent evolution; i.e., if the parent planetesimal was in 
a planet's resonance, is the particle still in resonance, 
and if so how are the parameters describing the reso- 
nant libration affected, which in turn tells us about the 
spatial distribution of such dust grains (e.g., W03). For 
the latter, the particles' subsequent evolution is not so 
important, since these orbits are hyperbolic. The more 
important issue is where those dust grains are most often 
created, and so where their hyperbolic orbits start, and 
that is an issue which depends on the collision rate of the 
parent planetesimals. The former bound grains are con- 
sidered in the remainder of this section, while the latter 
hyperbolic (blow-out) grains are considered in ^ 

2.1. Numerical Technique 

To derive the orbital parameters of small dust grains 
created in the destruction of planetesimals previously 
trapped in resonance by a migrating planet a three step 
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process was used: (i) first a population of parent plan- 
etesimals tliat were trapped in resonance by a migrating 
planet was defined; (ii) then the orbital elements of dust 
particles created in collisions between those planetesi- 
mals was worked out; (iii) then the dynamical evolution 
of those particles was followed to quantify the effect of ra- 
diation pressure on their resonance libration parameters. 
This process was repeated until the libration parameters 
could be determined for any size of dust grain associated 
with planetesimals trapped in a planet's 3:2 and 2:l(u) 
resonances. The new dust libration parameters were then 
used to work out the spatial distribution of these grains. 

The numerical technique employed in this paper is sim- 
ilar to that of W03. That is, numerical simulations were 
performed in which the dynamical evolution of 200 mass- 
less objects and 1 planet of mass Mpi was followed using 
the RADAU fifteenth order integrator program (Ever- 
hart 1985). All bodies are assumed to orbit a star of 
mass Mj, = 2.5Mq. For some runs the adjustment de- 
scribed in W03 was employed that results in a constant 
planet migration, Op/; this adjustment was modified so 
that the migration rate decreased at a linear rate to zero 
at the end of a migration. The "objects" were either 
planetesimals or dust particles for which their dynamical 
evolution is also affected by radiation pressure and the 
P-R drag force characterised by the parameter p. 

2.1.1. Parent Planetesimal Distribution 

To define a population of planetesimals trapped in res- 
onance, runs were performed in which a planet was made 
to migrate through a disk of 200 planetesimals at a rate 
starting at hpi but decreasing to after a time tmig- 
At the start of the integration the planetesimals had 
eccentricities, e, chosen randomly from the range to 
Gmax = 0.01, and were randomly distributed in semima- 
jor axis between a narrow range oi and 02. Their in- 
clinations, /, were chosen randomly from the range to 
emaa;/2 rad, and their arguments of periastron, d), longi- 
tudes of ascending node, fi, and longitudes. A, were each 
chosen randomly from the range to 360°. 

The planet had all of these angles set to zero at the 
start of the integration, with a semimajor axis apu. This 
meant that the planet ended up at a semimajor axis of 
o,pi2 = cLpii + a.pitmig/'^'- The planet's migration rate was 
chosen so that all of the planetesimals were trapped in 
the resonance being studied. Trapping probabilities are 
a strong function of migration rate, with slower migra- 
tions resulting in higher probabilities, and these were de- 
rived in W03. Since it was the strongest 3:2 and 2:l(u) 
resonances that were being studied in this paper, the 
migration rate was chosen so that it was as high as pos- 
sible while still trapping all of the planetesimals into the 
appropriate resonance; slower migrations could have re- 
sulted in trapping into higher order resonances. The re- 
lationship between apu and oi was determined by the re- 
quirement that the resonance being studied was in front 
of the planetesimals at the start of the integration; the 
location of the p -\- q : p resonance is given by: 

a, = ap,[(p + (z)/p]2/3, (4) 

although the finite width of the resonance also had to 
be taken into consideration. The width of the initial 
planetesimal distribution, 02 — 01, was set to ensure that 
the planetesimals did not encounter the resonance at the 
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same phase, thus biasing the distribution of resonant an- 
gles, (j) (W03). The extent of the migration is best de- 
scribed by its impact on the mean semimajor axes and 
eccentricities of the planetesimal population. Higher ec- 
centricities mean that more migration has taken place by 
an amount that can be determined from the eccentricity- 
semimajor axis relation = [q/ {p + q)]ln{a/ai) (W03). 

2.1.2. Initial Dust Distribution 

The initial dust distribution was taken directly from 
the orbits of the planetesimal population at the end of 
the run. In general each of these planetesimals would 
have a different collision rate and so some of the 200 
planetesimals could contribute more to the overall dust 
population, however this effect was not taken into ac- 
count and 200 dust particles were produced with exactly 
the same positions and velocities as the parent popula- 
tion. Different sets of runs were then performed for dust 
particles all of which have the same radiation pressure 
coefficient, (3. The initial orbital elements of the dust 
grains in each set of runs was determined from equations 
(PJ-(|S1). No further planet migration was assumed when 
considering the evolution of the dust grains' orbits. 

2.1.3. Dust Orbital Evolution 

Even though the parent planetesimals were in reso- 
nance with a = Or, the dust particles may no longer be 
in exact resonance for two reasons. First, the location of 
the resonance has changed because a dust particle moves 
slower than a planetesimal at the same semimajor axis 
because it sees a less massive star. Thus to get the same 
ratio of orbital periods a dust particle must be orbiting 
at: 

ard = (l-/3)'/'a,. (5) 

Second, the semimajor axis of the particle has changed 
as described in eq. Q Taking only terms to first order in 
eccentricity and /3, it is possible to show that the particle 
is at a semimajor axis which is offset from the resonance 
by a factor 

Aa = Od — Ord ~ 0^/3(4/3 ± 2e), (6) 

where the ±2e term indicates whether the particle was 
created when the planetesimal was at pericentre (/ — 0) 
or apocentre (/ — 180°) respectively. In other words, 
smaller particles end up further from resonance by a fac- 
tor oc /3, and particles released at pericentre also end up 
further from resonance than those released at apocentre. 
This does not necessarily mean that the particles are no 
longer in resonance, however, since resonances have finite 
width. 

A particle is said to be in resonance if its resonant 
argument (j) is librating rather than circulating, where 

= (p + q)^d - p\i - qi^d, (7) 

and libration can be characterised by a sinusoidal oscil- 
lation 

(j) = (j)^ + A(j> sin 2TTt/t^ (8) 

of period t^ and amplitude A(j), about a centre (j)m ■ Upon 
creation, the resonant argument of the dust particle is 
different to that of the parent, since both the pericen- 
tres of the orbits, and the longitudes of the particles 
within those orbits have changed (although in practise 
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the change in pericentre has the larger effect on the 
change in (jj). Overall the change in (f) is small unless 
a particle's /3 is close to, or larger than the eccentricity 
of the planetesimal e. Since e must be non-negligible 
to cause observable structure (W03), and particles with 
even moderate values of /3 are found to no longer be in 
resonance (see later), this effect is not considered further. 

It is also worth pointing out that the action of P-R 
drag is expected to cause the resonant argument of the 
dust particle to librate about a centre that is slightly off- 
set from 180° (or the appropriate centre for the case of 
the 2:1 resonance, W03), since this is required for the 
resonant forces to impart angular momentum to the par- 
ticle to prevent its inward orbital decay. The new centre 
can be derived by making the sum of the inward decay 
due to P-R drag (eq. 23 of Wyatt et al. 1999) and the 
semimajor axis variation due to resonant forces (eq. 14 
of W03) to zero, giving 

(j),n - 180° CX ^/fl, (9) 

where ip = Py^ M^^/a and fi = Mpi/Mi,. The parameter 
V' is equivalent to the parameter 9 in W03 in that it 
is the ratio of the dust particle's migration rate to its 
orbital velocity, which determines the angle at which the 
resonance is encountered.^ 

To ascertain the impact of radiation pressure on the 
particles' resonant arguments, the particles orbits were 
integrated for a sufficient amount of time for several li- 
bration periods to be completed, and so for the libration 
parameters to be fitted for each particle, and the mean 
parameters for each population of 200 dust grains to be 
determined, i.e., (0m}, (At/)), and (i^). Again the libra- 
tion periods of the parent planetesimals were derived in 
W03, and while the periods of the dust grains do differ 
from that of the parents, this provides a good enough 
idea to determine the required integration times. Since 
the oscillation also results in an oscillation of the par- 
ticles' semimajor axes, the amplitude of this oscillation 
was also determined, Aa, as was the mean for the popu- 
lation of planetesimals, (Aa). 

2.2. Results 

2.2.1. 3:2 Resonance 

The parameters used in the runs to characterise the 
libration parameters of dust grains created in the de- 
struction of planetesimals previously trapped in the 3:2 
resonance of a migrating planet are given in Tabled In 
this paper, a set of runs means those with the same star, 
planet and planetesimals parameters, but different values 
of /3 (corresponding to different sized dust grains). The 
first result is that for each set of runs, as /3 is increased, 
all of (0m), (A(/)), and (Aa) also increase, while (t^) de- 
creases (see Fig. 0. That is, for low (3 (large grains) the 
dust particles remain in resonance, but one with a libra- 
tion centre that is further from 180°, and with a higher 
libration width and range. All of the increases are linear 
with /3, although there is a turnover in — 180° for 
large (3. The runs found that particles remain in reso- 
nance as long as A0 < 180°. Larger particles (those with 

^ To put the angle i/i into perspective, a particle with /9 = 0.01 
at 1 AU from a IMq star would meet the orbital velocity at an 
incident angle of 0.4 arcscc (and would have -0 = 0.01). 



TABLE 1 

Parameters of the star (M*), planet (Mpi), 

PARENT planetesimals (a AND e), AND DUST 
PARTICLES (/3) FOR RUNS CHARACTERISING THE 
LIBRATION PARAMETERS OF DUST GRAINS CREATED IN 
THE DESTRUCTION OF PLANETESIMALS PREVIOUSLY 
TRAPPED IN THE 3:2 RESONANCE OF A MIGRATING 



PLANET. 




Mpi,M(s 


a, AU 


e 


/3 




2.5 


10 


78.6 


0.28 


0.0001 - 


0.004 


2.5 


10 


69.0 


0.20 


0.0001 - 


0.004 


2.5 


10 


64.4 


0.13 


0.0001 - 


0.004 


2.5 


30 


78.6 


0.28 


0.0001 - 


0.005 


2.5 


100 


78.6 


0.28 


0.0001 - 


0.007 


2.5 


300 


78.6 


0.28 


0.0001 ■ 


- 0.01 


0.5 


10 


78.6 


0.28 


0.0001 - 


0.006 


1.0 


10 


78.6 


0.28 


0.0001 - 


0.004 


1.5 


10 


78.6 


0.28 


0.0001 - 


0.004 


2.5 


10 


39.3 


0.28 


0.0001 - 


0.006 


2.5 


10 


196.5 


0.28 


0.0001 - 


0.004 



higher (3) are no longer in resonance and for such grains 
(f> circulates, i.e., undergoes a monotonic decrease. 

Different sets of runs were made varying each of the 
star, planet and planetesimal parameters in turn to show 
the way the libration of the dust grains was affected by 
each of these parameters (e.g., W03). This showed that 
the displacement of the libration centre from 180° de- 
pends on tp/fJ. as predicted in equation (j^. There is a 
turnover for high ijj/ fi, however the linear portion of the 
curve can be well fitted by: 

0m - 180° = 7740V'/Ai, (10) 

where n is in units of /Mq , and this line is shown on 
Fig- 13^ along with the results of all the runs. Since the 
libration width is also high for high '(/'/yLt, the significance 
or meaning of the turnover is not clear, and may be an 
artifact of the numerical method e.g., due to deviations 
from perfect sinusoidal oscillation in this regime. 

The increase in libration width could also be 
parametrised to explain the results of all runs: 

A0 = 16° + 88000/?M"°■^ (11) 

and this line is shown on Fig. along with the results 
of all the runs. In other words, the increased libration 
width is not dependent on radial distance of the planet 
from the star, or on the star's mass except in the ratio of 
the planet mass to the stellar mass, fi. The constant in 
this equation is indicative of the libration width inher- 
ent in the planetesimal population (W03) . The increased 
libration range could also be explained by the relation de- 
rived in equation © with the modification that higher 
libration ranges result from different planet masses ac- 
cording to: 

Aa/a,. = (4/3)/3 + 2 x 10"^, (12) 

and this result is represented in Fig.[5]3- A libration range 
above that given in equation © is expected if the parti- 
cle's new semimajor axis is not the peak in its libration. 
This is more likely to be the case for higher planet masses, 
since the libration range is higher for such planets (Mur- 
ray & Dermott 1999). 

The most important relationship in this paper is that 
given in equation Hll|) , since this can be used to estimate 
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Fig. 1. — Evolution of the resonant arguments, for different sized dust grains originating from the break-up of the same planetesimal in 
one of the runs for each resonance, and its implications for the spatial distribution of those grains. The top row is for the 3:2 resonance, for 
which dust grains exhibit a sinusoidal oscillation of <j) with a libration width that increases for smaller grains until the particles fall out of 
resonance and ij> circulates. The bottom row is for the 2:l(u) resonance, for which a similar progression is found, except that the libration 
is no longer sinusoidal (see text for discussion). The runs for both resonances correspond to parent planetesimals with eccentricities ~ 0.3 
that were trapped in resonance with a SOA/q planet which migrated 45-60 AU from a 2.5Mq star. The far left plots show the path of 
resonant orbits in the frame corotating with then mean motion of the planet at equal timesteps for an orbit with an eccentricity of 0.3; the 
resonant arguments, <p, determine the orientations of the loopy patterns of these orbits. 

the /3 above which particles fall out of resonance, which 
occurs when particles have > 180°. In other words, 
grains are still in resonance as long as (3 > Pcrit, where: 

Pcr^t - 2 X (13) 

Since this corresponds to large grains, the black body 
approximation can be used to estimate the grain size this 
corresponds to using the relation D = 0.4(_Li/Af*)//? in 
/im (assuming a density of ^ 2700 kg m^'^; Wyatt et al. 
1999) to give: 

D,„t = 200(L,/A/^)Ai-°-^ (14) 

in /im. For example, dust grains arising from the de- 
struction of the Plutino population in the Kuiper Belt 
(objects in 3:2 resonance with Neptune) remain in the 
resonance as long as they are larger than ~ 50 /im. 

Equations H10|l and Hll|) also specify the spatial distri- 
bution of dust grains of different sizes, since Fig. 6 of 
W03, a modified version of which is reproduced here in 
Fig. n showed how the angle cj) determines the orienta- 
tion of the loopy pattern that a resonant object's orbit 
makes in the frame rotating with the planet. For this 
paper the only important parameter describing the dis- 
tribution of the resonant arguments of the dust grains 
is the increased libration width, since one implication of 
equation l|13(l is that the offset of the libration centre 
from 180° is always small: even the smallest grains still 
in resonance have 4>m — 180° < 15° y/M^Jafj.^^-^ , so this 
can be ignored unless the planet is of very low mass and 
orbiting close to the star (such as dust grains trapped 
in resonance with the Earth, Dermott et al. 1994). The 
increased libration width of the orbits of resonant dust 
grains causes the clumps they make to be azimuthally 
smeared out. This is illustrated in Fig. |21 which shows 
the distributions in the frame rotating with the planet of 
dust grains of different sizes in one set of runs. 

Fig. O also illustrates how particles that are no longer 
in resonance have an axisymmetric distribution, and as 

^ To make these figures, the locations of the dust grains were 
coadded at all times in the frame rotating with the planet, and the 
integration was long enough to cover a large number of libration 
periods so that the structure was independent of the length of the 
integration (100,000 year integrations are shown in Fig.l^. 
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Fig. 2. — Libration parameters for dust particles in the runs 
of Table [Tl i.e., for dust originating from planetesimals trapped in 
the 3:2 resonance of a migrating planet: (a) displacement of the 
libration centre, ipm, from 180°; (b) libration wid t h, A (f>; and (c) 
libration range, Aa. The fits from equations 1101 - 1121 are shown 
with dotted lines. 
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Fig. 3. — Face-on view of the surface density distribution of dust grains of different sizes (characterised by the parameter /3) created in the 
collisional destruction of planetesimals trapped in resonance with a SOA/q planet which migrated 45-60 AU from a 2.5Mq star: (top) dust 
from planetesimals in the 3:2 resonance with eccentricities of ^ 0.28; and (bottom) dust from planetesimals in the 2:l(u) resonance with 
eccentricities of ~ 0.32. The planet's location is shown with both a white plus and a black diamond, and its orbital motion is anticlockwise; 
the fields of view of the far right plots (for /3 > 0.5) cover twice the area of the other (/3 < 0.5) plots. The colour table is scaled linearly to 
the maximum density in the density in each plot. 



expected this occurs for particles with (3 > (3crit (eq. • 
While the particles are not in resonance, the resonant ar- 
gument is still relevant to this discussion, since the pat- 
tern of the particle's orbit in the frame rotating with the 
planet is still close to that shown in Fig.Q] The particles' 
new semimajor axes mean that there is a monotonic de- 
crease in the value of (f), which if additional perturbations 
due to resonant forces are ignored can be estimated to 
be: 

(/.d = -720°p/3(l ± 1.5e)/V, (15) 

where p = 2 for the 3:2 resonance. This means that 
the loopy pattern rotates so that the the clump in front 
of the planet's motion approaches the planet (with the 
other clump 180° away). The timescale on which a non- 
resonant particle has its pericentre at the same longitude 
as the planet can be estimated from equation (|15(l . since 
for this to happen would have to decrease from 180° 
to 0° which occurs on a timescale of 



ts 



18O°/|0d| =0.25V/[p/3(l±1.5e)]. 



(16) 



This is an important timescale because if the parti- 
cle's eccentricity is high enough, a close approach to the 
planet is possible causing the particle to be scattered 
onto a more eccentric and inclined orbit. ^ Since af- 
ter scattering the particle's pericentre (or apocentre) re- 
mains close to the orbit of the planet further close ap- 
proaches and scattering can then ensue. This scatter- 
ing process was studied by following the simulations of 
(3 = 0.02 and 0.2 grains coming from the planetesimals 
in the 3:2 resonance shown in Fig.Olfor 350 Myr, though 
excluding grains once their semimajor axis dropped be- 
low 40 AU (either because such grains would be scattered 
by interior planets or evolve due to P-R drag on to the 
star). These simulations showed that some grains were 
inserted quickly into the resonances outside the 3:2 res- 
onance (some P = 0.02 grains were trapped in the 8:5 
resonance, while the P = 0.2 grains populated a large 
number of resonances including the 2:1, 5:2, 3:1, 7:2, 4:1), 

* Note that in the absence of close approaches, such as might be 
the case for dust coming from planetesimals with low eccentricity, 
it can be shown that P-R, drag would cause particles to reattain 
the resonance of their parent planetesimals on a timescale that is 
independent of particle S1Z6: tpr — 1070a2/M, m years. 



in which case they remained there until P-R drag forces 
had increased their eccentricities to the maximum value 
for the resonance at which point the libration width in- 
creased until the resonance was unstable, a process which 
took ~ 30 — 50 Myr. The non- resonant grains were scat- 
tered relatively quickly (on timescales as low as tgcat), 
but remained in the system to undergo further scatter- 
ing events. These grains were eventually excluded fol- 
lowing scattering events which put their semimajor axes 
below that of the planet, and once P-R drag had further 
reduced their semimajor axes below 40 AU, a process 
which took a total of ~ 10 — 30 Myr. 

The resulting spatial distribution of the grains is not 
only axisymmetric, but also more radially and vertically 
extended than that of the parent planetesimals. The 
resonances did introduce a small dumpiness in the re- 
sulting distribution, however due to the variety of res- 
onances populated by a small fraction of particles, this 
is expected to have a minimal effect on the distribution, 
resulting in only a slight underdensity at the location of 
the planet for these grains. The scattering timescale itself 
(eq. is relatively short, since grains with /? = 0.5 to 
Pcrit encounter the planet on timescales of tscat = 0.25ip/ 
to 62.5/i~°'^tp/. While this is likely to be significantly 
shorter than the collisional timescale, the simulations 
showed that the dynamical lifetime of the grains may be 
much longer than their collisional lifetime. The scatter- 
ing timescale given in eg. 1 161 should thus be considered as 
the timescale on which the particle distribution becomes 
axisymmetric. 

2.2.2. 2:1 Resonance 

The parameters used in the runs to characterise the 
libration parameters of dust grains created in the de- 
struction of planetesimals previously trapped in the 2:1 
resonance of a migrating planet are given in Table |21 For 
these runs, the migration was set so that planetesimals 
were trapped only into the 2:l(u) resonance (i.e., with 
librating about ^ 270° causing a clump of material ~ 90° 
behind the planet's motion). The results of the runs for 
the 2:1 resonance are broadly similar to those for the 3:2 
resonance in that for increasing /3 (decreasing grain size), 
the libration width increases leading to clump smearing 
and eventually an axisymmetric spatial distribution. 
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TABLE 2 

Parameters of the star (Af*), planet (Mpi), 

PARENT PLANETESIMALS [a AND e), AND DUST 
PARTICLES {/3) FOR RUNS CHARACTERISING THE 
LIBRATION PARAMETERS OF DUST GRAINS CREATED IN 

THE DESTRUCTION OF PLANETESIMALS PREVIOUSLY 
TRAPPED IN THE 2:l(u) RESONANCE OF A MIGRATING 
PLANET. 







a, AU 


e 


/3 


0.8 


10 


79.4 


0.44 


0.0001 - 0.1 


0.8 


30 


79.4 


0.44 


0.0001 - 0.1 


0.8 


100 


79.4 


0.44 


0.0001 - 0.1 


0.8 


300 


79.4 


0.44 


0.0001 - 0.1 


2.5 


30 


95.3 


0.32 


0.005 - 0.02 


2.5 


300 


95.3 


0.32 


0.01 - 0.03 


2.5 


300 


31.8 


0.32 


0.01 - 0.03 


0.25 


30 


95.3 


0.32 


0.01 - 0.03 


2.5 


300 


95.3 


0.43 


0.01 - 0.03 


2.5 


300 


95.3 


0.20 


0.01 - 0.03 


2.5 


300 


95.3 


0.12 


0.005 - 0.03 



However, one significant difference was tliat as tfie par- 
ticles /3 was increased, the libration of (j) could no longer 
be described as a sinusoidal oscillation. This is illustrated 
in Fig. n which shows the evolution of (j) for different 
sized dust grains coming from the same planetesimal in 
one of the runs, as well as Fig. O which shows the spatial 
distributions of different sized dust grains coming from 
the same poulation of planetesimals. As /3 is increased, 
first the libration becomes asymmetric in that more time 
is spent at low values of (/>; this means that the clump 
smears out asymmetrically about 270°, causing a spa- 
tial distribution weighted toward the anti-planet direc- 
tion. Increasing /3 further leads to the particle switching 
between performing half a libration in each of the 2:l(u) 
and 2:1(1) resonances. The distribution of such grains 
has three defining features: a gap at the location of the 
planet; concentrations just in front of and behind the 
planet (from the extreme points of the libration); and an 
overdensity in the anti-planet direction, caused by reso- 
nant forces which slow down the evolution of (/) at ~ 180° 
(this is just noticeable in the (3 = 0.01 run of Fig. but 
dominates the evolution of some grains). Finally the par- 
ticle is no longer in resonance at which point (j) circulates 
rather than librates resulting in an axisymmetric distri- 
bution. 

The non-sinusoidal oscillation makes characterisation 
of the structure more problematic than for the 3:2 reso- 
nance, and for the purposes of this paper the runs were 
used to derive the (3 for which the libration width is so 
high that grains with f3 > (3crit sue no longer in reso- 
nance. The result is that the stability of grains from 
planetesimals in the 2:1 (u) resonance is remarkably sim- 
ilar to those from the 3:2 resonance. In fact, (3crit can 
be well approximated by equation (|13|l as illustrated in 
Fig. 131 since the cut-off has the same dependence on AIpi, 
Mi, and a, and occurs at the same size of dust grain. One 
difference with the 2:l(u) resonance is that there is a de- 
pendence on particle eccentricity, in that (3crit oc e^-0±0-5.^ 
i.e., dust from higher eccentricity planetesimals remains 
in resonance down to smaller grains. Equation H13|l is 
thus valid for the 2:l(u) resonance when the eccentric- 
ities are ~ 0.3. However, since high eccentricities are 
required to cause significant structure (W03), the differ- 
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Fig. 4. — Collision rates amongst populations of planetesimals 
trapped into resonance with a 30Mg) planet that migrated 45-60 
AU from a 2.5A/*: (left) 3:2 resonance; (right) 2:l(u) resonance. 
The collision rate (bottom) is the product of the volume density 
of cross-sectional area of planetesimals (top) and the relative ve- 
locity of collisions (middle). Volume densities and collision rates 
have been normalised to the maximum for all planetsimals. Each 
point represents the result for one of 40,000 planetesimals in the 
disks, and is plotted against with the longitude of the planetesimal 
relative to the longitude of the planet. 

ence caused by this effect is likely to be small and equa- 
tion (|13|l is used in this paper as a good approximation 
for the cut-off. 

3. DISTRIBUTION OF BLOW-OUT GRAINS 

The distribution of grains that are blown out of the 
system by radiation pressure as soon as they are created 
depends on where those grains are most often created. In 
order to quantify the structure of the disk comprised of 
grains created in the break-up of large planetesimals, the 
following model was devised. As input, this model took 
the positions and velocities of a population of planetesi- 
mals that had previously been trapped into resonance by 
a migrating planet. For comparison with the simulations 
of fjSJ simulations were considered that placed planetesi- 
mals in the 3:2 and 2:l(u) resonances of a 3OM0 planet 
which migrated 45-60 AU from a 2.5Mq star. However, 
rather than using a full integration to determine the out- 
come for 200 planetesimals, these populations were simu- 
lated using the approach of W03, since this allowed large 
numbers of planetesimals to be simulated without re- 
quiring lengthy integrations (40,000 planetesimals were 
considered here) . Then for each planetesimal, the rate at 
which it collides with other planetesimals was determined 
using the fact that this collision rate is proportional to 
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Fig. 5. — Surface density distribution of /3 = 1 dust grains 
originating from the destruction of planetesimals trapped into res- 
onance with a 3OM9 planet that migrated 45-60 AU from a 2.5Mt: 
(left) 3:2 resonance; (right) 2:l(u) resonance. The radial distribu- 
tion of surface density is shown in the top plot, and the azimuthal 
distribution at 80, 400 and 2000 AU from the star is shown in the 
bottom plot. All distributions are scaled to peak at 1, except the 
planetesimal distributions which peak at 1/3. 

the product of the volume density of cross-sectional area 
of nearby planetesimals and the relative velocity of their 
colHsions (Opik 1951; Wyatt & Dent 2002). These two 
factors were derived by considering the number of plan- 
etesimals within a defined radius (4 AU in this case) 
of the planetesimal, as well as the mean of the veloci- 
ties of those planetesimals relative to the planetesimal in 
question. These factors, as well as the resulting collision 
rates, are plotted in Fig. ^ 

Fig. ^ shows that as well as a higher volume density in 
the clumps (as expected), there is also a higher relative 
velocity of collisions in the clumps. This means that the 
collisions are not only more destructive, but also occur 
at a higher rate (by more than an order of magnitude) 
than outside the clumps. Thus a planetesimal is most 
likely to be destroyed at the locations in its orbit when 
it is inside one of the clumps (i.e., near pericentre), and 
this is where the trajectories of the blow-out grains are 
most likely to start. This means that the distribution of 
these blow-out grains will not be axisymmetric. 

To determine the spatial distribution of the blow-out 
grains, the evolution of the planetesimal population was 
considered for a total time which depended how far from 
the star the dust distribution needed to be considered 
(6 planet orbits were considered in this case allowing 
the distribution to be determined out to at least 2000 
AU). Since the planetesimal distribution is constant in 
the frame rotating with the planet, the positions and 
velocities of the planetesimal population at subsequent 
timesteps were assumed to be the same as those of the 
planetesimal population at the end of the migration, but 
rotated with the planet; the collision rate of those plan- 
etesimals was as determined in Fig. 0] In each timestep, 
a fixed number of dust grains was introduced, with ini- 
tial positions and velocities that were the same (suitably 
rotated) as those of planetesimals chosen randomly ac- 



cording to the collision rates of the planetesimals. The 
evolution of those dust grains was then integrated for the 
remaining timesteps under the action of stellar gravity 
and radiation pressure. Equilibrium was achieved imme- 
diately by assuming the planetesimal that produced any 
given dust grain was equally likely to have created one 
a time equal to the total integration time ago; i.e., the 
integration "starts" with dust grains far from the star. 
The result is a movie which shows how the planetesimal 
and dust population evolves with time. Since the dust 
distribution is constant in the frame rotating with the 
planet, this was also coadded in the frame rotating with 
the planet to determine the spatial distribution of those 
grains. This is shown in the far right plots of Fig. O 
for 13 — 1 and 10. It is also shown more quantitatively 
in Fig. which shows histograms of the radial and az- 
imuthal distributions of the grains for /3 = 1. 

The radiation pressure blow-out grains form spiral 
structure which emanates from the clumps; i.e., for the 
3:2 resonance there are two spirals that start at longi- 
tudes ±90° from the planet, and the 2:1 resonance ex- 
hibits one spiral that starts ~ 90° behind the planet. 
Note that in these simulations dust grains were also cre- 
ated when the planetesimals were outside the clumps, 
just in much lower quantities. The spirals unwind in the 
direction opposite to the planet's motion. This arises 
because the pattern speed rotates with the planet (i.e., 
relatively fast), and while the dust grains can start with a 
comparable velocity to the planet's orbital motion (par- 
ticularly since they are created near pericentre), their 
azimuthal motion is soon much slower than that of the 
planet as they recede from the star. Fig. [S] shows how 
the radial distribution of the (3=1 grains falls off c>c 
outside the region where they are being produced, as ex- 
pected for grains on hyperbolic orbits; a similar fall off is 
seen for /3 = 10 grains for which the surface density drops 
oc r~^'^. This figure also shows how the spirals becomes 
less pronounced (i.e., the contrast between the density 
in and out of the spiral is lower) at large distances from 
the star as the spiral structure diffuses at large distances. 
Further it illustrates how the single spiral of the 2:1 reso- 
nance is maintained out to greater distances than that of 
the 3:2 resonance because it takes longer for this to merge 
with the nearest winding. In fact the spiral structure of 
blow-out grains will diffuse much faster than illustrated 
on Fig. [SJ since the tightness of the winding of the spiral 
is determined by the factor /?: dust grains with higher 
f3 are accelerated out the system faster than those with 
lower P meaning that the resulting spiral is less tightly 
wound. Representative simulations were performed that 
made the assumption that the dust grains produced had 
P chosen randomly from the range 1-10. Spiral structure 
close to the clumps is still well pronounced, and only 
becomes more diffuse at 400 AU. 

The distributions shown for /3 > 0.5 in Fig. are only 
valid for a disk comprised only of the planetesimals con- 
sidered in that simulation (i.e., all trapped in the same 
resonance with a narrow range of eccentricities). This is 
because the structure is determined by the collision rates 
of those planetesimals. A real disk is not just comprised 
of planetesimals in one resonance, and the interaction of 
the resonant plantesimals with planetesimals that are ei- 
ther non-resonant, or in other resonances, will affect the 
collision rates. To assess this, we also performed a rep- 
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resentative simulation of the distribution of dust grains 
expected as a result of the planetary migration described 
in W03 to explain the structure of the sub-mm emission 
from the Vega disk. In this simulation planetesimals were 
trapped into a variety of resonances (though mainly the 
3:2 and 2:l(u)), and many remained on non-resonant or- 
bits. The order of magnitude increase in collision rate 
is still apparent in the clumps, and so two armed spiral 
structure dominates the distribution of blow-out grains. 
Since the two clumps are of different magnitude it is also 
worth noting that one of the spirals is more dense by a 
factor greater than simply the ratio of the clump densi- 
ties, since the production rate of those particles is propor- 
tional to the product of the clump density squared and 
the relative velocity of the collisions (which is slightly 
higher in the brighter clump because of the superposi- 
tion of the 3:2 and 2:l(u) resonances there). 

The model discussed so far in this section is only valid 
for the P > 0.5 dust grains that are created by the de- 
struction of large planetesimals. Such dust grains could 
also be created in the destruction of intermediate sized 
dust grains that are still on bound orbits. As Fig. 13 
showed, those intermediate sized grains could have a dis- 
tribution that is axisymmetric. The blow-out grains aris- 
ing from the intermediate sized grains, with (3 > f3crit, 
would be expected to have a radial distribution similar 
to that shown in Fig. |S1 but with an axisymmetric az- 
imuthal distribution. 

4. DISCUSSION 

4.1. Summary of Grain Populations 

In §SJ21 and 121 the distribution of small dust grains re- 
sulting from the destruction of planetesimals that are 
trapped in resonance with a planet were studied show- 
ing how they differ significantly from that of the plan- 
etesimals themselves. Here the results are summarised. 
There are three distinct grain populations each of which 
exhibiting a different spatial distribution: 

• (I) Large grains with j3 > Pcrit (i-e., typically > 
a few mm) have the same clumpy resonant dis- 
tribution as the planetesimals, albeit one which is 
slightly azimuthally smeared out for (3 = (0.2 — 

f )/^crzi5 

• (II) Moderate sized grains with j3crit < /3 < 0.5 
(i.e., typically a few /xm to a few mm) are no longer 
in resonance and have an axisymmetric distribution 
which is also more radially extended and vertically 
broadened than that of population I grains. Such 
grains may also have a short lifetime due to the in- 
creased chance of a close encounter with the planet; 

• (III) Small grains with (3 > 0.5 (i.e., typically < 
a few /im) are blown out of the system by radia- 
tion pressure immediately on creation and so have a 
density distribution which falls off as r oc r~^, how- 
ever the structure of these grains can be further di- 
vided into two subclasses: (Ilia) grains produced 
in the destruction of population I grains that ex- 
hibit trailing spiral structure which emanates from 
the resonant clumps, and (Illb) grains produced 
from population II grains that have an axisymmet- 
ric distribution. 



4.2. Predictions for Multi-waveband Imaging 

Observations in different wavebands are sensitive to 
different sizes of dust grains. Thus a disk could exhibit 
different structures when imaged in different wavebands, 
providing the two wavelengths sample different grain 
populations. This has the potential to provide a valu- 
able observational test of models which explain clumpy 
structure seen in sub-mm observations of debris disks 
as due to resonant trapping of planetesimals. It is also 
relevant to ask whether the sub-mm observations would 
sample dust grains with a spatial distribution that is sim- 
ilar to that of their parent planetesimals (i.e., population 
I grains), as was assumed in W03. 

The size of grains sampled in different wavebands de- 
pends to some extent on the composition of the grains, 
but is most strongly dependent on the size distribution 
of grains in the disk. At present this size distribution is 
not predicted by the models. The most simple assump- 
tion is that, since the dust is produced in the coUisional 
destruction of planetsimals, the size distribution is the 
same as that expected in an infinite coUisional cascade 
with n{D) oc D~^'^ down to the radiation pressure blow- 
out limit (Wyatt & Dent 2002). Although in practise it 
has been shown that an abrupt cut-off to the size distri- 
bution at the blow-out limit would cause a wave in this 
distribution (Thebault, Augereau & Beust 2003). This 
simple assumption can be used to explore the relative im- 
portance of grains in populations I and II, but does not 
make any predictions about the quantity or observabil- 
ity of population III grains. It also neglects the fact that 
the dynamical lifetimes of population II grains may be 
lower than their coUisional lifetimes due to the increased 
chance of scattering by the planet, and their coUisional 
lifetimes may be affected by the broadening of their spa- 
tial distribution. 

However, this simple size distribution has been shown 
to provide a good fit to the spectral energy distributions 
of the emission from several debris disks (Wyatt & Dent 
2002, Sheret, Dent & Wyatt 2004). The resuUs pre- 
sented in Fig. 5 of Wyatt & Dent (2002) can also be 
used to work out which grain populations we would be 
seeing in the Fomalhaut disk for different wavebands (as- 
suming of course that some of the planetesimals in this 
disk are trapped in resonance with a planet, which has 
not been proved). Since just 5% of the sub-mm emission 
comes from grains either < 300 /zm or > 20 cm in diam- 
eter, with roughly equal weightings in this range toward 
different bins of log particle diameter, this implies that 
sub-mm observations would indeed be dominated by the 
population I grains that have a similar spatial distribu- 
tion to that of the planetesimal population. However, 
this also implies that, since less than 5% of the 25-100 
fim emission comes from grains that are bigger than 6 
mm with weightings that favor the lower end of the size 
range, such observations would be dominated by popu- 
lation II grains and would have an axisymmetric spatial 
distribution. 

Thus for disks where this size distribution holds, and in 
which clumps are seen in sub-mm images, a model which 
interpreted the clumps as the result of planet migration 
and the consequent trapping of planetesimals into the 
planet's resonances could be tested by seeing if mid- to 
far-IR imaging of the disk shows it to be axisymmet- 
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ric. By consideration of equation H14() it is possible to 
infer that the wavelength at which the transition from 
clumpy to smooth structure occurs is indicative of both 
the mass of the perturbing planet, with the transition 
shifted to shorter wavelengths for more massive planets, 
and the spectral type of the star, with shorter wavelength 
transitions for lower mass stars. In other words multi- 
wavelength imaging could also reveal the mass of the 
planet. In fact Fomalhaut's disk appears asymmetric at 
a range of wavelengths (Holland et al. 2003; Stapelfeldt 
et al. 2004; Marsh et al. 2005; Kalas et al. 2005), favor- 
ing a mechanism that affects grains of all sizes such as 
an offset centre of symmetry due to the non-circularity 
of a perturbing planet's orbit (Wyatt et al. 1999; Marsh 
et al. 2005; Kalas et al. 2005). 

To assess the importance of the population III grains 
relative to that of the populations I and II grains de- 
scribed above, a more detailed model of the size distri- 
bution of the disk would be required which is beyond the 
scope of this paper, although some considerations are 
discussed in 



4.3. Implications for Vega 

Vega's dust disk has been observed to be clumpy when 
imaged in the sub- mm (Holland et al. 1998), yet more 
recent observations show its structure to be axisymmet- 
ric at mid- to far-IR wavelengths (Su et al. 2005). This 
does not in itself validate the model of W03, since Su 
et al. also showed that the disk is much more extended 
at these wavelengths than in the sub-mm. They devised 
a three component model to explain the radial emission 
distribution observed from 25-850 /im which is comprised 
of grains that are 4, 36 and 430 fini in diameter. The 
two smallest sizes have an optical depth distribution that 
falls off T oc and are small enough that (3 > 0.2. 
Such grains dominate the mid- to far-IR images and are 
interpreted as grains that are in the process of radia- 
tion pressure blow-out (i.e., population HI grains), while 
the larger grain population, with a density distribution 
peaked around 100 AU is required to fit the 850 /im im- 
ages (and must be predominantly population I grains to 
produce the observed clumpy structure). Su et al. in- 
terpreted this as evidence of a recent collision having 
occured in the disk, since the blow-out grains are short- 
lived and it is unfeasible for the observed outflow to have 
been continually replenished over Vega's 350 Myr life- 
time. However, they did not attempt to explain the az- 
imuthal structure in the sub-mm images, or lack thereof 
in the short wavelength structure. Here the model of Su 
et al. is expanded to turn the 3 component model into 
one with a size distribution extending across all sizes. 
The aim is to provide a more realistic description of the 
disk which can be used to question whether, in the ab- 
sence of a model that predicts the size distribution in 
the W03 model, the relative quantities of populations I, 
II and HI grains are physically plausible in the context of 
that model given our understanding of size distributions 
in coUisional cascades, or if another model needs to be 
sought to explain the clumpy sub-mm structure. 

The size distribution in Vega's disk is assumed to be 
defined by different power laws in three size ranges con- 
sidered to represent populations I, II and III grains. Su 
et al. already showed that the distribution of population 
HI grains follows a size distribution of n(D) cx D^^ ''^ in 
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Fig. 6. — Model for the size distribution of dust in the Vega 
disk: (a) Cross-sectional area as a function of particle size, with 
boundaries indicated between the different populations I, II and 
III. For population III the area is the total out to 1000 AU; (b) 
CoUisional lifetime of different sized grains in the disk, assuming 
the coUisional properties of weak ice from Wyatt & Dent (2002), 
and an average collision velocity of 1.4 km/s and a conversion of 
cross-sectional area to volume density that assumes this area is 
spread evenly around a torus extending 86-200 AU from the star for 
populations I and II, and only considering the area of population III 
grains in the same region; (c) Contribution of different grain sizes 
to the total Hux in the different Spitzer and SCUBA wavebands; 
(d) Spectral energy distribution of emission from the Vega disk 
model (solid line) showing the contribution of different populations 
to that spectrum (dotted, dashed and dash-dot lines). Fluxes are 
excess fluxes after photospheric subtraction from the IRAS Faint 
Source Catalog, Spitzer (Su et al. 2005), and SCUBA (Holland et 
al. 2005). The level of the photospheric emission is shown with 
the thick solid line. 

the range '^4 — 50 /im and to have a distribution t oc 
for r > 86 AU. This distribution is maintained here, but 
the minimum and maximum grain sizes, as well as the 
total cross-sectional area, are left as variables. The pop- 
ulation I and II grains are fixed to lie at 100 AU, the peak 
in the optical depth distribution found by W03, and the 
size distributions of both populations are assumed to fol- 
low that of a coUisional cascade with n{D) oc D~^-^, but 
with different amounts of cross-sectional area in each of 
the populations. The division between populations I and 
II was set at 1 mm, a changeover size consistent with a 
planet of size 40 - 50M® (eq. Following Su et al., 
astronomical silicates were assumed for the particles' op- 
tical properties. Since Su et al. already showed that this 
fits the surface brightness distribution, it remains to fit 
the spectral energy distribution (SED), where the emis- 
sion from population HI grains is taken to be that out 
to 1000 AU. An additional constraint was set to deter- 
mine the relative contributions of the population I and II 
grains, which is that both of these populations combined 
contribute ~ 50% of the total fiux at 160 /im as found in 
the modeling of Su et al. 

The size distribution which fits the observed emission 
spectrum (and by analogy with the modeling of Su et 
al., also the surface brightness distributions), is shown 
in Fig. 1^. The contribution of different size grains in 
the distribution to the fluxes in the Spitzer and SCUBA 
wavebands is shown in Fig. Et, and the modeled emis- 
sion spectrum, showing the contribution of the different 
populations, is shown in Fig. Because of the way 
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the emission effieieney falls off for population II grains at 
longer wavelengths, the condition of equal contributions 
from bound and unbound grains to the 160 /im flux (and 
the need to fit the 850 /im flux) results in a size distribu- 
tion with a level of population II grains some 3-4 times 
lower than that expected in an idealised coUisional cas- 
cade (e.g., that assumed in Wyatt & Dent 2002). This 
also results in the population II grains contributing just 
~ 25% of the 850 /im flux, in line with the observa- 
tion that this image is clumpy (Holland et al. 1998). 
Increasing the relative amount of population 11 grains 
both increases their contribution to the 850 /im image 
and of bound grains to the total 160 /im flux. It is this 
result, along with the quantity of material in population 
III relative to populations 1 and 11, that needs to be ex- 
plained within the context of coUisional cascades. The 
results for the boundaries defining population 111 grains 
are less important, since these values depend on the op- 
tical properties of the grains which are not considered 
here. ^ 

To answer whether the distribution of Fig. |^ is a re- 
alistic, or even expected, size distribution in Vega's col- 
lisional cascade, Fig. 03 shows the lifetimes of the dif- 
ferent grains due to collisions. This plot exhibits some 
well known features: starting at the largest sizes, the 
coUisional lifetime of grains is reduced as smaller sizes 
are approached due to the greater quantities of cross- 
sectional area of smaller objects in a coUisional cascade. 
The lower numbers of population II grains, however, 
cause an increase in the lifetimes of the population I 
grains that should have been destroyed by this popu- 
lation (Thebault et al. 2003). The population II grains 
have a lower lifetime than the population I grains because 
these are destroyed in collisions with the large quantities 
of blow-out grains in this distribution (Krivov, Mann & 
Krivova 2000). However, the lifetime of the population 
HI grains themselves due to collisions is much lower than 
their blow-out time indicating that only a small fraction 
of the population HI grains are themselves destroyed in 
collisions on their way out. This means that it may be 
possible to explain the break in the size distribution be- 
tween populations I and II, as inferred from the SED 
modeling, by the destruction of population II grains in 
collisions with those in the process of blow-out by radi- 
ation pressure (Krivov et al. 2000), and the knock-on 
effect induced by the lower levels of such grains in the 
distribution at larger sizes (Thebault et al. 2003). This 
effect may be further accentuated by a lifetime for popu- 
lation II grains that is even shorter than their coUisional 
lifetimes of ~ 100,000 years which may be caused by 
interaction with the planet which occurs on 1000 year 
timescales fea. I16|l . 

However, this still leaves the question of why there are 
so many population HI grains. This modeling finds a 
similar mass loss rate to Su et al. for this population of 

2M0/Myr, which, as Su et al. point out, implies that 
the cascade must have been initiated relatively recently. 

^ The change-over from population II to III in this model occurs 
at 60 fim, which is at /3 = 0.17. A cut-off at such a low /3 may 
be real, since a cut-off at /3 < 0.5 is expected if, as shown in 
Fig.|l] most grains are created at the pericentre of highly eccentric 
orbits (since such grains are put onto hyperbolic orbits provided 
/3 > (1 — e)/2). However, this small discrepancy could equally be 
removed by changing the optical properties used in this model. 



However, the lack of evidence for non-axisymmetry in 
the far-IR emission at large distances from the star seen 
by Su et al. poses problems for the W03 model. This is 
because the mass of population II grains is comparable 
to that of population III grains at ^ 2 x lO^^M®, yet 
the population II grains are only destroyed in collisions 
on timescales of ^ 100,000 years whereas population III 
grains are removed on 1000 year timescales. This means 
that, unless there is some mechanism that turns popu- 
lation II grains into population HI grains on 1000 year 
timescales, then at most 1% of the population III grains 
can be of type Illb. The remainder must be population 
Ilia grains, and the collision rates indicate that this is 
not unreasonable: for the distribution assumed in Fig. El 
the mass in population I grains (0.6M^) is processed in 
collisions at a rate ~ 7Af0/Myr. While most of this mass 
is likely redistributed within population I rather than lost 
to population III, it is noted that this rate would be much 
larger if the distribution had been assumed to extend to 
sizes larger than 10m. Thus it is possible that the re- 
quired loss rate could be achieved with just a few % of 
the mass of population I objects being put into blow-out 
grains in destructive collisions as long as the population I 
distribution extends to large enough objects. This means 
that the mid- to far-IR emission should exhibit spiral 
structure rooted in the clumps seen in the sub-mni. Lim- 
its on the circularity of the emission detected by Spitzer 
were not discussed in detail in Su et al. , but it is possible 
that these images were not of sufficiently high resolution 
and/or sensitivity (or are too confused by the point-like 
photospheric emission) to rule out the presence of pop- 
ulation Ilia grains. If this is the case, then the model 
predicts that at high resolutions the disk should exhibit 
spiral structure when imaged at far-IR and mid-IR wave- 
bands, e.g., when imaged at 25 /^m using a coronagraph 
with MIRI on the JWST (Wright et al. 2003). A more 
detailed confrontation of the model with the Spitzer and 
SCUBA observations of the Vega disk is left for a future 
paper (Wyatt, Su, Rieke, Holland, in prep.). 

5. CONCLUSIONS 

This paper shows how the distribution of small dust 
grains resulting from the destruction of planetesimals 
that are trapped in resonance with a planet differ from 
that of the planetesimals themselves, both in terms of 
their orbital characteristics and consequently their spa- 
tial distributions. Three different grain populations are 
identified based on grain size: population I grains that 
are large enough to remain in the resonance of the parent 
object (and so have a clumpy distribution); population II 
grains that, due to radiation pressure, are still on bound 
orbits, but are no longer in resonance (and so have an ax- 
isymmetric distribution); and population HI grains that 
are removed from the system by radiation pressure on 
short timescales (and so have a distribution that falls off 
oc 7'^^). Subclasses are defined for population 111 grains 
based on the population designation of the parent ob- 
ject: population Ilia grains originate in the destruction 
of population I grains and exhibit trailing spiral struc- 
ture emanating from the clumps; population lllb grains 
originate in the destruction of population II grains and 
have an axisymmetric distribution. 

The fact that a planetesimal belt is made up of parti- 
cles from all populations and subclasses, each of which 
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has its own (quite different) dynamical and spatial distri- 
bution implies that observations in different wavebands 
can be dominated by different populations and so ex- 
hibit different morphologies. Adoption of a simple col- 
lisional cascade size distribution with no blow-out (pop- 
ulation III) grains implies that if some planetesimals in 
the disk are trapped in resonance with a planet then sub- 
mm observations would trace the distribution of those 
planetesimals through population I grains (validating the 
approach of Wyatt 2003 in modeling the SCUBA obser- 
vations of Vega), but that mid- to far-IR observations 
would trace population II grains and so have an axisym- 
metric distribution. The wavelength at which the transi- 
tion occurs from clumpy to smooth structure is indicative 
of the mass of the planet. 

The size distribution of Vega is modeled in the light 
of recent Spitzer observations which show that signifi- 
cant quantities of population III grains are present (Su 
et al. 2005). This shows that there is a significant lack 
of population II grains in this distribution (a factor 3-4 
under that expected in a collisional cascade). Analysis of 
the collisional lifetimes indicates that this may be due to 
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the destruction of these grains by those that are in the 
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argued that unless there is some mechanism that is de- 
stroying the population II grains on 1000 year timescales, 
then the population III grains, and so the mid- to far-IR 
images of the Vega disk, should exhibit spiral structure 
emanating from the clumps seen in the sub-mm images. 
Detection of such structure may be possible with Spitzer 
or with MIRI on the JWST, and would confirm the in- 
terpretation of the morphology of Vega's disk in terms 
of planetesimals trapped in resonance with a planet or- 
biting at 65 AU (Wyatt 2003) as well as indicate the 
direction of its motion. 

Multi-wavelength imaging thus provides a method of 
confirming models interpreting clumps in debris disks as 
indicative of planetesimals trapped in resonance with an 
unseen planet. Such images can also provide information 
on the mass and direction of motion of the perturbing 
planet, and do not require the decade timespans of multi- 
epoch methods for confirming these models by checking 
for orbital motion of the clumps (e.g., Ozernoy et al. 
2000). 
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